Draft version February 1, 2008 

Preprint typeset using LAT^^X style cmulateapj v. 14/09/00 



(N 
O 

o 

(N 



A PERSPECTIVE ON THE CMB ACOUSTIC PEAK 
T.A. Marriage^'^ 

Draft version February 1, 2008 

ABSTRACT 

CMB angular spectrum measurements suggest a flat universe. This paper clarifies the relation between 
geometry and the spherical harmonic index of the first acoustic peak {£peak)- Numerical and analytic 
calculations show that ipeak is approximately a function of Qk/^m where fix and Qm arc the curvature 
{^Ik > implies an open geometry) and mass density today in units of critical density. Assuming 



Hk/^m <C 1, one obtains 



-peak 



llV3 



where and 



are the scale factor 



at decoupling and radiation-matter equality. The derivation of tpeak gives another perspective on the 
widely-recognized ^m-^k degeneracy in flat models. This formula for near-flat cosmogonies together 
with current angular spectrum data yields familiar parameter constraints. 

Subject headings: cosmic microwave background — cosmological parameters 
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1. INTRODUCTION 

It has long been recognized that the first peak in the 
CMB angular spectrum provides information about the 
curvature of the universe (Doroshkevich, Zel'dovich, & 
Syunyaev 1978; Bond et al. 1994; Kamionkowski, 
Spergel, & Sugiyama 1994; Efstathiou & Bond 1999; 
Cornish 2000; Weinberg 2000). The data are now in. 
TOCO, BOOMERanC, MAXIMA, and DASI have mea- 
sured the peak position (Miller et al. 1999; Netterfield 
et al. 2001; Lee et al. 2001; Halverson et al. 2001). 
Analyses of the data strongly suggest a flat universe (Hu 
etal. 2001; Jaffeet ah 2000; Stompor et al. 2001; Pryke 
et al. 2001; Dodelson & Knox 2000). The MAP satellite, 
cosmic variance limited through £ « 600, should make the 
definitive measurement in the near future (Page 2000). 

As new data resolve higher order peaks, attention shifts 
to the physics at angular scales beyond that of the first 
maximum. The object of the present analysis is to clar- 
ify the physics derived from the position of the first peak, 
hereafter called "the peak index" or Ipeak ■ The size of the 
sound horizon r^* at an angular diameter distance Da* 
to decoupling determines the peak index. This is widely 
recognized and serves as a starting point. In Sections 2, 3, 
and 4, numerical and analytic calculations yield the peak 
index as a function of CIm and Qk- All results are checked 
with CMBFAST (Seljak & Zaldarriaga 1996). In Section 
5, a simple formula for ipeak, applicable to low Qk/^m 
universes, is developed alongside geometric and classical 
interpretations of the formal results. It is found that the 
peak index approximates a function of Qk/^m- Although 
our analysis grounds itself in the familiar concepts of Da* 
and Tg*, the results are not widely recognized. While the 
physical effects responsible for ipeak are understood, the 
interplay that gives the flx/^M dependence is not intu- 
itive and holds a number of surprises. We end the inves- 
tigation by using the functional form of ipeak and current 
angular spectrum data to obtain parameter constraints re- 



sembling those of more sophisticated treatments. 

Throughout this work, the {^Imj^k) plane serves as the 
parameter space. Other quantities enter, though the re- 
lation between ipeak and CIk/^m is relatively indepen- 
dent of them. To be concrete, we take ilBh"^ = 0.02, 
consistent with nucleosynthesis (Buries, NoUett, & Turner 
2001). The redshift of decoupling, z*, is taken as 1400 
from the Saha equation. The Hubble constant assumes 72 
km/s/Mpc in accord with the HST Key Project results 
(Frcedman et al. 2001). 



2. ANGULAR DIAMETER DISTANCE 

The comoving Friedmann-Robertson- Walker 
metric for a closed universe (Qk < 0) may be 

as 



(FRW) 
written 

ds^ = dr/^ ~ {dx^ + sin^ ^rfS^) (1) 

where 77 denotes conformal time, = —l/flx, and dS^ 
is the line element of a unit two-sphere. All spacetime 
intervals are given in units of Hq^. The angular diam- 
eter distance Da to an object of proper length ds at co- 
moving distance Rx {drj = dx = 0) is defined so that 

dS =1 ds I /Da. 

Da = Rsmx{a) (2) 
is the solution to the Friedmann equation: 

-da2 



where 
dx'- 



(3) 



Flat and open geometries are treated analogously. 

Numerical Result Da* can be evaluated numerically. 
The distance to redshift 1400 is computed for models 
across the {Q,m,^k) plane using cosmography routines 
from David Hogg (Hogg 1999). 

Analytic Result To complement the numerical result. 
Da* can be estimated analytically by setting JIa, Oeg, a* to 
zero in equation (3), and then using Mattig's solution (cf. 
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Numerical Result Analytic Result 




Fig. 1. — The angle subtended by the sound horizon at decoupling. AGs = fs^/Da*- Data favors AO^ = 0.60° (Knox, Christensen, & 
Skordis 2001). In both numerical and analytic plots, equation (7) gives the sound horizon Ts*- In the numerical plot, the angular diameter 
distance to last scatter Da* is computed using the cosmography code (Hogg 1999). In the analytic result, the classical equations for Da* 
are integrated with Q.f^,aaq,a* = 0, and the resulting formulae are used without assuming CIm + = 1 (see eq. [4]). AGg from either 
method approximates a function of Q,k/^m- The bending of the contours in the analytic result arises from radiation-dominated dynamics 
before last scatter. In the numerical result, explicit inclusion of in the Priedmann equation balances the radiation effect: AGg contours 
are straight lines converging on zero in the full calculation. At far right is illustrated the comoving coordinates of a closed FRW universe with 
time and azimuthal angular coordinates suppressed. Inspection of the drawing shows AGs = rs*/Rsvax* where the comoving sound horizon 
at decoupling Ts* and the comoving distance to decompling Rx* are calculated using Priedmann dynamics. 



(Peebles 1993)). 



Da* 



x5(7) 



where 7 = 2 | fix \ /^m, and S=l for a flat universe. If 
is non-zero, then S takes the form 



5(7) = 7172^; 

(3 = 7+l + \/(7Tl) 



VI -(7-1)^ 



1. 



nK>o 



(5) 



3. THE SOUND HORIZON 

Sound travels through a tightly coupled photon-baryon 
system with speed Cg 

Cs = , ^ (6) 
n/3(1 + Q) 

where Q is ZpE/^P-y, and ps and p^ are the energy den- 
sities of baryons and radiation, respectively. The sound 
horizon at decoupling rg* in comoving coordinates is then 



= / Cs—da 
Ja=o da 



^/l + Q* + ^Q* Qeg 



(Hu & Sugiyama 1996). Note that curvature does not af- 
fect sound dynamics before decoupling. In any geometry, 
Ta* will be the same. 

4. THE HORIZON ANGLE AND THE PEAK INDEX 

Da* and Vg* give the angular size of the horizon: 

^5* 



Ae« = 



intersect the origin of the (f^M, ^k) plane. The angle sub- 
tended by the sound horizon, and so the position of the first 
(4) peak, approximates a function of Qk /^M ■ A well known 
corollary to this general statement is that a peak corre- 
sponding to fix = should be insensitive to variations in 
(Bond et al. 1994; Hu et al. 2001). 
To check whether this simple analysis agrees with the 
standard model, we compare the above results to those 
fromCMBFAST (Seljak & Zaldarriaga 1996). CMBFAST 
calculates the CMB angular spectrum."' Then, Ipeak from 
the spectrum multiplies A8s from Figure 1 to give the 
constant of proportionality a — ipeak x AO^. As shown 
in Figure 2, the numerically derived a increases from 110 
{flM = 0.14) to 125 (Qm = 0.74) and has weak if any 
dependence on curvature. The near constant graphs of 
a suggests a well-defined agreement between the peak in- 
dices from CMBFAST and those derived from this work's 
numerical and analytic calculations (see middle frame of 
Figure 2). 

5. DISCUSSION 

To gain intuition about AG^, consider the closed comov- 
ing universe (eq. [1]) illustrated in Figure 1, where time 
and azimuthal angular coordinates have been supressed to 
produce the familiar two-sphere geometry. In this picture, 
CMB photons follow great circles, and, assuming a small 
angle, AQg derives from inspection: AOg = rs*/Rsmx*, 
which is exactly equation (8). 

To better understand the parameter dependence of rs* , 
take the sound speed Cg (eq. [6]) to be constant at 9/10^/3. 
(With Qsh^ =0.02 and radiation density derived from the 
COBE FIRAS measurement, the sound speed decreases at 
a near constant rate (dcg/da « constant) from c/\/3 at 
/g\ a=0 to four-fifths that value at decoupling (Fixen et al. 
1997).) The sound horizon is then 
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{y/a* + aeq - ^Ja^ (9) 



Values for AGg corresponding to the numerical and ana- 
lytic results for Da* (see Section 2) are plotted in Figure 1. y 
Curves of constant AOg approximate straight lines, which 

4 CMBFAST inputs are [Q-b = 0.04, f2„ = 0,i?o = 72,Tc^i, = 2.726, Y^e = 0.24, N^{massless) = 3.04, N,,{massive) = O,recfast,no reion, 
scalar only,primordial index l,adiabatic). 



Cg / —da = — p== 
Jo da 5V3S2m 
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Fig. 2. — The relation a = AGs x ^peak and parameter constraints. The peak index from CMBFAST times AGs of numerical and 
analytic calculations gives a (left frame). The curves for a are averaged and reinserted into the equation ipgak = AGg/a to generate the 
corresponding graphs of peak index (middle frame). Contours of (peak iii the {Qm,^k) plane are calculated using equations (9) and (10) 
(right frame). To coincide with the values of AGs and tpeak derived from current data, the contour plot assumes a = 125(Ri 0.6° x 210) 
throughout the plane (Knox, Christensen, & Skordis 2001; Knox & Page 2000; Hu et al. 2001). The dashed contours show constraints from 
current estimates of tpeak E'ld from the projected MAP satellite measurement. 



where the final equality follows from the Priedmann equa- 
tion (3) with fix = f^A = 0. 

To obtain a simple formula for ipeak, expand Da* to first 
order in 7 = 2 | Qk \ /^m' 



2 + 



K 



(10) 



The expansion applies to models with a low curvaturc-to- 
mattcr ratio. Coincidentally, such models are favorc;(l by 
experiment, so that the resulting formula for ip^ak is use- 
ful. Given equations (9) and (10) and a value of a from 
Figure 2, 



aDa 



^peak 



IIV3 



This near-flat approximation is plotted in Figure 2 where 
tteq « 2.4 X lQ~^{VtMh'^)~^ ■ For flat models, the distance 

to last scatter (eq. [10]) scales as ^j^' . One may intuit 
that self-gravitation leads to smaller cosmological separa- 
tions. This same "gravitational" effect decreases the dis- 
tance between the big bang and last scatter, and therefore 

— 1/2 

Tg* (eq. [9]) also scales with an overall factor of fl;^ ' . In 

— 1/2 

equation (11), the dependence of Da* cancels that 

of rs*. This cancellation helps explain the Qm-^a degen- 
eracy of flat cosmogonies. Furthermore, the Qk/^m term 
in equation (11) is proportional to the curvature times the 
area between light rays in a two-dimensional representa- 
tion of an Qk ~ universe (e.g. Figure 1). This suggests 
that one may think of {^m,^k) dependence of ipeak in 
near-flat spacetimes as resulting from light curving like 
the geodesies of a two-dimensional space of constant cur- 
vature. Finally, the effect of radiation on early universe 
dynamics and, in particular, on rg* is made explicit by the 
appearance of a,.g in equation (9). Radiation-dominated 
cosmological growth per expansion scale (drj/da) is less 



than that of matter- dominated dynamics. In low Om uni- 
verses, radiation brings last scattering even closer to the 
big bang and so shortens Vs* . Thus, the effect of radiation 
in the early universe is to spoil the pure / depen- 
dence of ipeak as manifest by the curved contours in the 
plot of equation (11) in Figure 2. The straightness of the 
f-peak contours is restored in the numerical result shown in 
Figure 1. This suggests that explicit inclusion of in the 
computation of Da* balances the Ugq dependence of r^* . 

6. CONCLUSION 

The peak index has long been recognized as an indicator 

of geometry. It is hoped that the present analysis sheds 
new light on ipeak- The peak index does not determine 
the magnitude of curvature, but rather the ratio of cur- 
vature to matter. A measurement of the peak's angular 
scale gives the precise geometry only if Six ~ 0, otherwise 
ipeak is a fmiction of Qk/^m- Furthermore, in deriving 
the I^m dependence of (.peak , unexpected cosmological 
cancellings were discovered. Particularly useful is the bal- 
ance of overall matter dependencies in r^* and which 
helps account for the flM-^A degeneracy in flat models. 
At the same time, however, it is remarkable that the ad- 
mittedly simple arguments of this work yield such a de- 
cisive cosmological indicator. Within the next few years, 
NASA's MAP satellite data should give Ipeak to cosmic- 
variance levels. This measurement will burn a sharp line 
of possible worlds across the {flM,^K) plane. 

This letter was begun at Cambridge as part of DAMTP's 
tripos. TM thanks Ofer Lahav and Daniel Wesley for early 
guidance, David Hogg for cosmography computer code, 
and Lyman Page and James Peebles for suggestions re- 
garding the final draft. TM is an NSF Graduate Research 
Fellow and is supported by NSF grant PHY-0099493. 
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